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O . Abstract 

(N ■ ^ . 

• For general relativistic equilibrium stellar models (stationary axisymmetric asymp- 

Q . totically flat and convection-free) with differential rotation, it is shown that for a 

C/^ I wide class of rotation laws the distribution of angular velocity of the fluid has a 

^sO ' sign, say "positive", and then both the dragging rate and the angular momentum 

. density are positive. In addition, the "mean value" (with respect to an intrinsic 

I density) of the dragging rate is shown to be less than the mean value of the fluid 

' angular velocity (in full general, without having to restrict the rotation law, nor the 

■ uniformity in sign of the fluid angular velocity); this inequality yields the positivity 

! and an upper bound of the total rotational energy. 

OS • 

O : PACS numbers: 04.40.Dg, 97.10.Kc, 02.30. Jr 
O ■ 

1j : I. INTRODUCTION 
u '. 



One of the most interesting of the relativistic effects produced by the rotation of a star 
is the dragging of inertial frames (also called Lense-Thirring effect). 1^ This has been 
classically described in terms of its local effects on gyroscopes and particles. However, 
^ I its cumulative effects on the motion of particles give a much simpler description: a lo- 
I cally non-rotating test particle, that is dragged along in the gravitational field of the 
star, has an angular velocity, as seen from a non-rotating observer at spatial infinity, 
which is named angular velocity of cumulative dragging, or rate of rotational dragging, 
shortly called dragging rate."^ ^ It is physically expected that, for isolated rotating stars 
in thermodynamic equilibrium, this dragging rate A has the same sign (rotation sense) 
as the fluid angular velocity fl, if this one has a "uniform" sign throughout the fluid (in 
the general differentially rotating case). Indeed, LindblonP and, independently, Hansen 
and Winicour (1977)^ seem to establish this result, however, without explicitly fulfill- 
ing the corresponding requirements when applying the called Hopf theorem (a maximum 
principle) to an elliptic operator in a certain domain, concerning the boundedness of its 
coefficients on the boundary of the domain, specifically on the axis of rotation, and the 
(and not C^) regularity of the metric functions across the surface of the star. 
Also, (assuming in the description above that the test particle does not collide with 
the star's matter if it goes through the star) one is tempted to conjecture, in principle 
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in the rigidly rotating case with f2 = const. = fi* > 0, that the dragging rate is bounded 
above by the fluid angular velocity, A < Q^,. And Hansen and Winicour (1975^ offer some 
proof of this (although with the same objection as above). 

In the general differentially rotating case, however, an analogous relation should not 
be expected to be so simple; firstly, because different portions of the star's interior could 
have opposite rotational motion about the same axis (assuming a convection-free fiuid), 
and, secondly (even if the fiuid angular velocity has a sign), because, due to the inte- 
grability condition of the equation of motion, the distribution of fiuid angular velocity, 
f2-profile, cannot be freely prescribed; instead, it is derived (together with the potential 
functions, integrating the field equations) once an appropriate rotation law is given. Most 
of the literature concerning numerical works on differentially rotating neutron stars make 
generally the ansatz for a certain rotation law which yields A < Q. Nevertheless, for a 
more general law such a relation is not so obvious. Hansen and Winicour (1977)P have 
made some attempts to give a result, however they needed the unphysical assumption 
that the star's matter occupies the whole space. 

One of the aims of this work is precisely to find general and physically reasonable 
assumptions on the rotation law of a differentially rotating stellar model, so that the 
dragging rate is (at each interior point) less than the fluid angular velocity, and, hence, 
the angular momentum density is positive (vanishing on the axis) provided the weak en- 
ergy condition is satisfled. For that matter we consider the time-angle fleld equation's 
component, which is elliptic and linear in the dragging potential A in coordinates adapted 
to the symmetries. The approach with the metric in these coordinates is attractive because 
the fleld equations become semilinear elliptic. Specially, they reduce to four (coupled) el- 
liptic equations for the four metric functions. One has however to control the coordinate 
singularities of the equations on the axis of rotation, but these can be treated mathemati- 
cally using the axial symmetry of the physical problem.- So handled, the elliptic equation 
in A writes in a "regular" form and has bounded coefficients; this allows us to apply a 
maximum principle to several differential inequalities, which, using the C^-matching on 
the star's surface and the asymptotic flatness of the metric, will lead to the mentioned 
and other interesting inequalities. 

More generally, as was conjectured by Thorne^ (p. 245), the mean value of A is ex- 
pected to be less than the mean value of Q. However, to my knowledge, there is in the 
literature no explicit and so general result in this direction. The other purpose of the 
present work is then to derive a "general" inequality on "mean values" with respect to 
a density function. In addition, related to this question is the concept of total rotational 
energy of the star. Hartl^has given bounds for this rotational energy in the slow rotation 
limit, which we aim here to generalize. 

The paper is organized as follows. After reviewing in Sec. II the model for a relativis- 
tic star which is rotating differentially. Sec. HI is devoted to handle the concerned fleld 
equation, elliptic and linear in the dragging potential, with special attention to the regu- 
larity and boundedness properties of the involved functions, as a preparation allowing us 
to apply the maximum principles (reviewed in the Appendixes) in Sec. IV, where inequal- 
ities concerning mainly the positivity of the dragging rate and of the angular momentum 
density are derived. In Sec. V a general "mean values inequality" is derived in full gen- 
eral; and the positivity and upper bound of the total rotational energy is established in 



2 



the general differentially rotating case. Finally, in Sec. VI the relevant results are brieffy 
summarized. 



II. MODEL FOR A DIFFERENTIALLY ROTATING RELATIVISTIC STAR 

The spacetime of an isolated rotating star in thermodynamic equilibrium within general 
relativity theory is generally represented by an asymptotically flat stationary axisym- 
metric 4-dimensional Lorentzian manifold with metric g = Qapdx'^dx^ satisfying 

Einstein's equations, 

Gap '■= Raf3 — -^RgaP = 871 Taf3 , (1) 

for the energy-momentum tensor of a perfect fluid, T^fs = {e +p)uaUf3 +pgai3, with 
4-velocity m", energy density e, and pressure p. Signature of g is here considered to be 
( — h ++). Since the star is isolated, the matter (perfect fluid) is confined in a compact 
region in the space (interior), with vacuum, Taf3 = 0, on the outside. 

We denote the two (commuting) global time and axial Killing vector fields® by ^ = 5^ 
and r] = d^, respectively, where = t labels the space-like hypersurfaces which are 
invariant under time translations, and x^ = (f) is the axial-angle coordinate around the 
axis of rotation, given hj rj = 0; {t,(f)) G Mx [0,27r[. The metric components will then 
only depend on the two remaining spatial coordinates, gap = gai3{x^,x^)- 

We shall assume that the fluid motion is purely azimuthal (non-convective) , i.e. the 
fluid 4-velocity is contained in the 2-surface spanned by the two Killing fields, (as 1-forms) 

u A ^ A r7 = (circularity condition) . (2) 

In that case it can be seerP that the 2-surface elements orthogonal to the 2-dimensional 
group orbits of the Killing fields are surface forming (the same holds in the vacuum region); 
and, consequently, the metric may be written in a form which is explicitly symmetric under 
the change (t, 0) (— t, — 0). In the 2-surfaces orthogonal to the orbits we can always 
introduce isotropic coordinates {x^,x^) = {p,z) without loss of generality, so that the 
metric can always be reduced to the standard fornP^ 

g = gapdx'^dx^ = -e^^dt" + e"'^ [p'e^''{d<p - A dtf + e^'^idp' + dz")] , (3) 

where the metric functions U, B, and A depend only on the (p, 2;)-coordinates of the 
"meridian plane". Here p and z are cylindrical coordinates at the asymptotically flat 
infinity, and, using the remaining freedom of conformal transformations in the meridian 
plane, we choose these coordinates such that p = represents the axis of rotation and 
(p, z) e Mq X M (denoting := {x G M | x > 0}). The metric functions pe^, U, and A 
can be written as invariant combinations of the Killing fields in the form 



p2e2B = _ det(((?^,)^,,=t,^) = -g(|, $,) g{r), -q) + g(^, r}f 

g(^7,'^) (4) 
g(^,^) ' 
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and they can be interpreted physically as follows: p e represents a sort of distance from 
the rotation axis (and, hence, B is, to some extent, a measure how far is p from being 
that distance); [/ is a generalization of the gravitational potential; and A is the angular 
velocity of cumulative dragging, or dragging rate. The remaining metric function is K, the 
conformal factor in the meridian plane. 

Throughout the following we shall denote the closure and the boundary of a set X by 
X and dX, respectively. We fix the notions 



/ = interior of the star := {(p, ;z) G R^^ x R | p{p, z)>0} cR^xR 
E = exterior of the star := (M+ x M) \ 7 c M([ x M (5) 
S = star's surface := 7 n^= 9/ C Mc^x R, 

/ and E open in the induced topology in M.q x R c R^; that means, although part of the 
axis (p = 0) is in / (and part in E), the only points of the axis which are in dl = S (and 
in dE = S U {oo}) are the poles, if they exist. The set / C Rq x R is supposed to be 
bounded and connected. Concerning the regularity of 5* = dl, we assume it satisfies an 
exterior sphere condition everywhere (cf. Definition in Appendix A). 

Within our star model, the matching conditions (from the interior and the exterior 
solutions) require that the pressure vanishes identically on the star's surface, p = on S. 
In the exterior {Tajj = 0) we have e = p = 0. Furthermore, e and p satisfy a barotropic 
equation of state in the interior, 

e = e{p) in / . (6) 

We assume the pressure p to be continuous with respect to the coordinates, and also 
p e{p) a continuous function, 

p G C°(R+ X R) , p e{p) G C°(R+) , (7) 

satisfying the weak energy condition^ 

e + p>0 (in R+ X R) . (8) 

(Notice, by the definition of the interior, if e > in J, as it is generally assumed, we 
shall have even e + p > in I, and, hence, condition (jHJ follows. In addition, since the 
equation of state is defined only in the interior, requirement (|7j) does not guarantee 
the continuity of e across the star's surface (where p = 0), namely, if e{p = 0) > 0, then 
a jump discontinuity of e across the star's surface occurs.) 

From the circularity condition Q on the fluid 4- velocity (in /), this is of the form 

d(j) 

u = u*(^ + flrj), where Q = — = — 

w- dt 

is the angular velocity of the fluid measured by a distant observer in an asymptotically 
fiat spacetime, and the fact that the 4- velocity u is a unit time-like vector field determines 
the normalization factor u*, such that g(u, u) = —1, i.e. 

(«*)-2 = e2^-pV(^-^)(^-^)' =: N, (9) 
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from wliere N = (m*) ^ > in /. Indeed, we do not allow that the velocity of light is 
approached somewhere, and, hence, even 

N > const. > in 7. (10) 

We consider a star rotating differentially with a distribution of angular velocity {rotation 
profile) Q = fl{p,z), a continuously differentiable function, 

n G c\7) . (11) 

However, the r2-profile of the fluid cannot be freely chosen, this shows up in the following. 
The integrability conditions of the field equations ((H), that is, the equation of hydro- 
static equilibrium T"^.^ = (from G"^.^ = 0) (where ';' denotes covariant derivative), 
particularly, its part orthogonal to the fluid 4-velocity u, reduces to the Euler equation, 

dp=~{£ + p)sL, (12) 

where a is the 4-acceleration of the fluid, a = VuU. Specifically, 

8i = dV + u^u^dn, V=^\nN, (13) 

u^u^ = p^e^'^^~^\VL — A)N^^. But the integrability condition of Eq. ()12|) taking into 
account (jH)) is da = 0; following then, from (fT3|l . diu^u^) f\ dVL = The special case 
where = const, is called rigid rotation (or uniform rotation). In general we shall have 
Q 7^ const., following then, 

M% = j^in) , (14) 

for some function JF, rotation law. By specifying the function J-'{Q) a specific model 
of differential rotation is obtained. (Note, since in the Newtonian limit m'm^ — ^ Q, 
Eq. fll4|l expresses the general relativistic generalization of the Newtonian "rotation on 
cylinders" theorem, Q = Gip"^) ). 

Further requirements on our stellar model are: 

a. the metric functions are (at least) two times continuously differentiable in the in- 
terior and in the exterior of the star, and continuously differentiable everywhere 
(cf. Note in Subsec. III.B), 

K, U,B,A e C\I) n C\E) n C\R+ X R) ■ (15) 

b. in order that the metric functions are symmetric with respect to the 2;-axis (p = 0) 
( "axisymmetric solutions"), and, hence, the metric (jHl), defined on Ai excluding 
the axis, can be extended to an at least axisymmetric tensor field in the whole 
spacetime J\4, we assume that 

as p ^ 0, dpK, dpU, dpB, dpA , (16) 

and, for completeness, also dpE, dpp ; 

c. finally, by the asymptotic flatness requirement, denoting V := {dp,dz), 

as := (p2 + z2)i/2 ^ oo, K,U,B,A-^0 and VK, VU, VB, VA ^ . (17) 

Notice, from regularity, in (fT5|l . and asymptotic flatness, (fT7j) . it follows, in partic- 
ular, that the metric functions and their derivatives are bounded,'^ 

|f/|, |A| < oo and \\VK\\, \\VU\\, \\VB\\, \\VA\\ < oo in R+ x R . (18) 
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III. THE ELLIPTIC EQUATION FOR THE DRAGGING RATE A 

A. The time-angle field equation component 

The {t(j)) component of Einstein's equation in these coordinates takes the fornP 

dppA + d,,A +-dpA+ {3VB - AVU , VA) = -^j^ ■ {Q - A) , (19) 

with ip'^ := 16n—{e + p) (> 0, by © and ^) , (20) 

where (.,.) denotes the Euchdean scalar product. Since, from condition (fTI)jl . v{p,z) = 
v{—p, z) for V = K,U, B, A, i.e. we are considering only axisymmetric solutions of the 
field equations, and since only "axisymmetric operations" appear in these equations, we 
consider the following transformation (in the spirit of Ref. 6) in order to avoid the coordi- 
nate singularity (of Eq. ()19|l ) on the axis of symmetry (2;-axis, i.e. p = 0). To this end we 
use the 5- lift of each function v = v{p, z) (on R^), for the metric functions v = K,U, B, A 
and also for v = Q,e,p, where the n-lift of v : R.q x M ^ R on flat R", axisymmetric 
around the defined as follows 

V ^ V such that v{x) = v{xi, . . . ,Xn) := v(^p = [x^^ + ■ ■ ■ + x^^iY^'^ , z = Xn) , (21) 

and, for every function v ■.W"' —>■ R, the meridional cut (in direction xi) ofv, 

V y-^ V such that t'(p, 2;) := ^(p, 0, . . . , 0, 2;) . (22) 

For axisymmetric functions, both operations are isomorphisms and inverse to each other; 
but the relevant properties of n-lift and meridional cut are that (a) they leave the regular- 
ity conditions and the norms invariant, (6) they commute with "axisymmetric operations" , 
in particular, with all operations in Eq. ()19p . like multiplication and scalar product, yield- 
ing specially (for n = 5) 

{Vv, Vw) = {Vv, Vw) , denoting V := (9i, . . . , ^5) {di = d^^, dij = d^^d^^ ) , 

and, remarkably, (c) they transform the operator dpp + dzz + ^^^^p ^ 2) into the fiat 
ra-dimensional Laplacian, and vice versa; having for n = 5 

3 ^ 
dppV + dzzV + -dpV = ^ diiV =: Av . 

P i=i 

Hence, with the 5-lift, Eq. (I19|l writes in the form 

AA + {3VB - AVU , VA) = -ip^ ■ {n - A) , (23) 
ijj'^ defined like ip"^, pOjl . but with 5- lifted functions (on R^). 

B. Regularity and boundedness of the metric functions 

Let us see how conditions (fT^ - (fTHj) transmit through the 5-lift. First, from conditions 
(US)) and (Uni) it follows 

k,U,B,A e C^(R^), (24) 
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because (for v = K,U,B,A) v e C\R+x R) and dpV ^ as p imply d G C\M.^). 

Note: In fact, as seen in Ref. 6, with the use of these mathematical tools (n-lift and 
meridional cut, for different numbers n), the elliptic system of field equations (P) may be 
regarded as a set of Poisson-like equations, where the nonlinearities (quadratic terms in 
the first derivatives of the metric functions) are contained in the inhomogeneous terms 
on the right hand side. Making the weak requirement that the metric functions and their 
derivatives are essentially bounded, v, Vv G L°°, since also e, p E L°° (by condition ((Zj) 
and e = p = in the exterior), and fl G (by (fTT|) ). we have that the right hand side 
is essentially bounded. Then, by the regularity of Poisson's integral,'^ (at least) v G C^'" 
for some a < 1; in particular, v E C^, i.e. fl^ . This justifies requirements (fT3|) and fTBj) 
in Sec. II. 

Combining (fT3j) and (fT^ we obtain also that the 5-lifted metric functions are class 
in the interior and in the exterior of the star (in MP). That is, denoting 

T := {(xi, . . . , X5) G I ((xi^ + . . . + x^Y^ , 0:5) G /} C (25) 
and, analogously, S and S, from E and S (cf. ©), respectively, we have for v = K,U, B, A 

V G c\i) n c\S) n c^W') . (26) 

The asymptotic flatness condition (fTTj) implies, through the 5-lift, that 

as R= \\x\\ = {xi^ + . . . + x^^y/^ oo, v^O and -> ; (27) 

but V G C\R^), that is, d G C^{R^) and Vv G [^^(MS)]^ , yielding, together with condi- 
tions pTj) . their respective boundedness, 

\v\ < oo and ||V{;|| < oo in . (28) 
C. Notation convention and roundup 

We have seen in Subsec. III. A that Eq. fl^ is equivalent to Eq. flTUj) through the 5-lift 
and the meridional cut, ()21|) and (j^^ for n = 5. Furthermore, the 5-lift leaves regularity 
and boundedness properties invariant; see Subsec. III.B. 

Convention: For simplicity in the notation, we omit throughout the following the symbols 
' ~ ' for all 5-lifted functions we use. (Once it has been seen how regularity and boundedness 
properties transmit from the functions defined on Rq x R to the lifted ones (on R^), and 
since they are equivalent in terms of positivity, and no explicit reference to the first ones 
will appear throughout the following section, this notation convention seems appropriate.) 

Accordingly, we write Eq. (PHj) in the form 
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in -A), 


with Lq a 


= AA + {3VB - AVU , VA) , 






:= 


167T—{e + p) >0 (=Oin^), 


and 


:= e - p e ''(^7 - A) > const. > m J, 



(29) 
(30) 
(31) 



where K,U,B,A : 



^ — > M, axisymmetric around the xs-axis. (Notice, Eq. ()29|1 is so 
defined in the whole spacetime, interior (fiuid) and exterior (vacuum), 2US = M^, but in 
the exterior ip"^ = {e = p = 0) and the vacuum field equation is recovered, Lq A = in 
£.) Also, we have ()26|) - (PH|) . i.e. (with the notation convention) 



K, U, B, A E c\i) n C^{£) n C 



1 /toSn 



(32) 



oo, K,U, B, A^O, 

VK, Vf/, V5, VA 



0. 



(33) 
(34) 



\K\, \U\, \B\, \A\ < oo 
and llVi^ll, ||V[/||, ||V5||, ||VA|| < oo in 



(35) 
(36) 



Equation (jSHl), i-e. 

writes then 



LA:= LoA-i/j^ ■ A = -ijj^ -n, 



LA = aij{x)dijA + bi{x) diA + c(x) A = g{x) 



(37) 



with = const. = 5ij (= 1 if z = j, and = otherwise) , 

bi = 3d,B-AdiU (Vi, j e {1, . . . , 5}), and 

C = (<Q)^ 

g = cn, 



(where repeated indices indicate summation from 1 to 5). The fiat 5-dimensional Lapla- 
cian A, in (jSOI), (ajj = ^ij), and hence L, is obviously strictly and uniformly elliptic 
everywhere. The coefficients hi are measurable and bounded functions everywhere, be- 
cause B and U are C\ (jH^ . and have bounded derivatives, (jHSl). On the other hand, 
for the coefficient c (cf. ^^), since (i) the metric functions are continuous, (jH^ . and 
bounded, (jH^j) : (ii) p is continuous everywhere, ((Tj), and has compact support; (iii) e is 
continuous in the (closed) interior X, from ((Zj); (iv) f2 is in particular continuous (in X), 
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(fTTj). and, hence, measurable; and (v) > const. > (also in X), (fTUj) : it follows that 
c = —tp"^ is measurable and bounded in the interior X, and, since i/j"^ = [e = p = 0) in 
the exterior S, and the boundary (the star's surface) dl = 5 is a set of measure zero, we 
have that the coefficient c is measurable and bounded everywhere. This will allow us in 
the following section to apply maximum principles in the classical and in the generalized 
sense to the operator L (and Lq); see Appendixes A and B. 

IV. BOUNDS ON THE DRAGGING RATE 

A. Positivity of the dragging rate 
Proposition 1 

// the distribution of angular velocity of the fluid is non-negative (and non-trivial) , then 
the dragging rate is positive everywhere. 

n>o, n^o =^ A>o. 

Proof. Consider the domain G defined by a ball in MP centered at the origin x = and 
of arbitrarily large radius a, 

G := B„{Q) C . (39) 

Since A is continuously differentiable in M^, cf. (jH^ . so is in particular in G; but A and 
continuous in implies that they are 2-integrable (are in L^) in G; consequently, 

A e iy^'2(G) nC^(G). (40) 

Hence, the strictly elliptic linear partial differential equation (in A) with measurable 
and bounded coefficients, pTjl . is satisfied in a generalized sense in G; see Appendix B. 
Remarkably, whenever ^2 > 0, Eq. ()37|1 yields the differential inequality 

LA<0 inG, (41) 

i.e. ^4 is a generalized supersolution relative to the operator L and the domain G. We pay 
now special attention to the behavior of A on the boundary: since the radius of the ball 
G, cr, is arbitrary, we can make it sufficiently large (a — >■ oo) such that, by the asymptotic 
flatness condition on A (cf. fl33|) ). A is arbitrarily small on 9G, 

lim A\qq = 0. (42) 

We first observe that A ^ const, (because, by (jl^ and A G G°(M^), would he A = 
const. = 0, which yields, by Eq. (jTfjl . f2 = 0, and we are assuming ^ 0). Hence, by the 
strong minimum principle, Theorem|3]in Appendix B, applied to the differential inequality 
(PT|) . A cannot attain a non-positive minimum at an interior point of G; using we 
conclude then A > in G, i.e. everywhere. □ 

Remark 1. A result analogous to Proposition Q holds with the opposite sense of the 
rotation; that is, if i7 < (fi ^ 0), then A < 0. This follows because Eq. fl29|l is invariant 
with respect to the simultaneous change of sign (f2. A) — > (— f2, — A). 
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B. Upper bound Q. Positivity of the angular momentum density 

Hereafter we discuss the sign of the difference Q — A. Remarkably, this determines 
the sign of u^u^, which, with assumption (jH)), is the sign of the angular momentum 
density, integrand of the total angular momentum, given by the "volume" integraP^ 
J = Jj 27rT^* {—gf^'^dx, where g = det(g) and = {e + p) u^u^. 

1. In the rigidly rotating case 
Proposition 2 

In the particular case of rigid rotation, with Vt = const. =: > 0, 

< A<n^ 

holds everywhere. As a consequence, in this case, u^u^, and, hence, the angular momentum 
density, is non-negative. 

Proof. We consider Eq. for f2 = const. = > 0, i.e. LqA = —ip'^ ■ — A), which, 
since the differential operator Lq, ()30|). is free from linear term, can be rewritten in the 
form 

LoiA - n,) = -i? ■ (a - A) , 
or, denoting again the differential operator L := Lq — ip"^ and defining 

w{x) := A{x) - (43) 

in the whole spacetime, x & X U S = M.^ (as already 5-lifted function; cf. Subsec. III. A), 

Lw = L{A - ^],) = Lo{A - fi,) - ip^ ■ {A - n^) = in . 

We have then the strictly elliptic linear (in w) equation 

Lw = (in particular) in G = B^{0) C , (44) 

where the radius a is arbitrary, with w G W^'^{G) n C^(G'), by (jHI and (jlHl). On the 
other hand, by the condition of asymptotic flatness on A (in ()33|)). A is arbitrarily small 
on dG, provided that a is sufficiently large, i.e. (jl^ : consequently, 

lim w\gQ = -fi* < . (45) 

cr — 5-00 

Since w ^ const, (because, by (jl5|l and continuity, would he w = const. = — ^2* in G, 
that is, A = const. = in G, which is not allowed, by Eq. (jHTjl . since here = ^7* > 0), 
applying the strong maximum principle. Theorem^ in Appendix B, to Eq. (@H), we get 
that w cannot attain a non-negative maximum at an interior point of G; hence, using 
w < in G (everywhere), i.e. A < everywhere. Moreover, A > everywhere, by 
Proposition H This establishes the conclusion of the proposition. □ 

Observe, in the static case, = 0, we would have Lw = and limg-^oo Mac ~ 0' 
following, by the strong maximum and minimum principles, w = 0, i.e. A = = 0; as 
expected, A = 0. 
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Remark 2. Likewise, if = const. = fi^ < 0, then > A > ^l^ everywhere, and, hence, 
the angular momentum density is non-positive. We obtain this by applying Proposition |21 
to the function A := —A, solution of Eq. fl^i for := —Q^ > (cf. Remark Q). More 
explicitly, the angular momentum density of a rigidly rotating stellar model has the same 
sign as the angular velocity of the fluid. Also, as a result, we have for a fluid rotating 
rigidly with Q = const. = Q^, 

< |A| < . 



2. In the general (differentially rotating) case 

In the following we shall assume that a function JF (to be specified) has been given, 
and we have a solution of the problem, that is, (four) metric functions, K, U, B, and A, 
and a fluid angular velocity distribution, Q, satisfying the (four) field equations (H)) (in 
particular, the elliptic equation for A, Eq. and Eq. (HH), u^u^ = J-'{fl). (Notice, 

in the interior, where the matter terms do not vanish {p > 0, e > 0), substituting into 
the equation of motion (fT2|l (with (fT3|l ) its integrability condition, i.e. Eq. (fT4|l . and the 
equation of state, Eq. (jH)), we obtain the pressure, p, and the energy density, e, as functions 
of p,U,B,A, and fi.) 

Remarkably, u^u^ may be written 



pV(^-^)(l]-A) 
g\n-A) 



■.^{g,n-A) with e:= pe^-^^ , (46) 



i-g-'in- Ay 



where, from (fTO|l . 1 — g'^{Q — AY = Ne > const. > in T. With the defined function 
(gni), Eq. (mi) writes 

^{Q,n- A) =J^{n). (47) 

Lemma 

Assume 

i. the function T : M — > M is strictly decreasing, and 

ii. 3 a constant Qc (|^c| < oo) such that J-'{Qc) = 0, 

then, at each interior point (in I), where Eq. (jTTj) . $ = JF, zs satisfied, the following holds 

A<n ^ A<n<n^ iA<n^) 
A>n ^ A>n>n^ > 

A = Q ^ A = Q = Q^ {A = Q,). (48) 

Note 1: Due to (i), fi^ (defined in (ii)) is unique. Also, observe, VL^ exists and coincides 
with the (constant) value of VL on the rotation axis, provided that part of the axis, ^? = 0, 
is in the interior, X, (i.e. if the rotating fluid does not have toroidal topology). This 
is because at points in {f? = 0} fl X 7^ 0, since $|g=o ~ ^ ^"^^ $ = in X, we have 
^{^)\q=o = 0) ^'^^ ^ is, by requirement (i), invertible; yielding fi|^=o ~ const. = ilc- 
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Note 2: Observe, if J" G and J" ' < 0, then, since > 0, Eq. ()47|) can be solved 
for fi, by virtue of the implicit function theorem, yielding Vt = Q{p, U, B, A); and, by the 
regularity of the metric functions, (jHSj), it follows in particular E (X) , requirement 

Note 3: It should be stressed that, since $ is an increasing function in Q, choosing 
the function JF strictly decreasing (requirement (i)), Eq. ()47|1 has a unique solution in Q 
("curve" solution with g variable). Specially, this makes likely the existence of functions 
fl, K, U, B, and A, solutions of the field equations and Eq. (fT^ . Indeed, in numerical 
works concerning differential rotation the ansatz for the jF-law J-'{Q) = Rq{Qc~^), where 
i?o is a free parameter describing the length of scale over which Q changes, is generally 
used, and they claim they have a solution. (See e.g. Refs. 15-17.) 

Proof. We consider a point x El where the metric functions and the fluid angular velocity 
are solution, in particular, with reference to Eq. (j47jl . the functions $ and valued at 
this point "intersect" each other, i.e. 

^{q{x) , fi(x) - A{x)) = J^{n{x)) (Vx E I) . 

From requirements (i) and (ii), it follows (at each interior point) sign(jF) = sign(i7c — ^)- 
As regards $ (at the interior point), on the axis {q = 0) it obviously vanishes, cf. 
following, from the relation $ = JF (in X), $ = JF = and, thus, Q = Qc on the axis. 
Outside the axis {g ^ 0) we have sign($) = sign(fi — A), and, hence, by $ = JF (in X), 
sign(Oc — f^) = sign(f2 — A) outside the axis. This yields (pSjl . □ 
We are now in a position to get one of the main results of this work in the general 
differentially rotating case, namely. 

Proposition 3 

// the T-lavj (in Eq. ()14|1 ) specifying the model of differential rotation is chosen such 
that 

i. JF : R — > R strictly decreasing, 

ii. 3 a constant flc (|^c| < oo) such that JF(fic) = 0, and 

iii. Qc > 0, 

then 

0<A<n<n^ znT; (49) 

in particular, u^u^ > 0, and, hence, the angular momentum density is non-negative. 
Moreover, 

< y4 < max < fi^ in £ . (50) 

Note: As remarked above, if the interior (fluid) contains points of the axis, then condition 
(ii) is already guaranteed, and VL^ is the constant value of Vt on the axis; cf. Note 1 in the 
previous lemma. See also Notes 2 and 3. And observe, requirement (iii) is in principle 
much weaker than f2 > 0, but, as seen in the conclusion of this proposition, > already 
follows. Furthermore, the fact that VL <VLc'va.X shows that in differentially rotating stars 
the core may rotate faster than the envelope, so that the core can be supported by rapid 
rotation before mass shedding is reached at the equator.'^ 

Proof. We divide the proof in four steps. 
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1st step: Let us see first A < Q in X. 

Suppose (to get a contradiction) A{xo) > fi(xo) for some Xq G X. We have seen, in 
the previous lemma, cf. (jlSj) . that this is equivalent to A{xo) > fl{xQ) > flc', and, 
hence, using hypothesis (iii), A{xo) > fic > 0. Therefore, by the continuity of A 
(indeed A G C^, cf. (IH^ ) and the asymptotic flatness (hmn^^n^oo ^ = 0, cf. (jHSl)), 
we infer that 

3 an open and connected neighborhood of xq, Mxq C M^, such that 

A>VL^ in A/;^ (51) 
and A = VLc on dNx^ ■ 

We distinguish two cases: 

case 1: Af^^nS = 0, that is, the neighborhood is contained in the interior, A4„ C X. 
Thus we have, again using the previous lemma, A > Q > flc > mMx^; particularly, 
il- A <0 in UxQ, and, therefore, by Eq. (jSHl), 

LoA>0 in7V;.(,CX. 

From dSl, in particular, A G ^^(X) n C°(X), and, hence, A G ^^(TV;,,) n C°(Ar^,,). 
And applying the weak maximum principle. Theorem ^ in Appendix A, to the 
operator Lq (on A) in Mx^-, we obtain that the maximum of A is reached on the 
boundary, i.e. 

max A = max A , 

contradicting (|HT|l . 

case 2: A^o H £^ 7^ 0. (Notice, here is included the case where xq G dl = S.) We 
denote 

Ix, = A^oHX CX 

and r = A4o n 5 . 

Observe, T 7^ 0, because A^o is connected, Xq G A^q fl X 7^ 0, and, by assumption 
in this case, A4o fl £^ 7^ 0. Moreover, T = Ixo ^ ^xq = dTx^ H dSx^. 

Thus, we have 

in the interior, from (jS^, 

A>n, in ix, u r c Ko (r c 9x,j 
A = n, on 9x,„ \ r c 9A4o ; 

and, applying the weak maximum principle (Theorem Q in Appendix A) to the 
differential inequality (cf. (jSH))) 

LoA>0 in X^o C X 
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(again using (gHl), A > > fi^ in J^^J, with A G ^^(J^J n C^iXxo), it follows 

max A = max A = max A : 

in the exterior, we have analogously, from ()5ip . 

A > fi, in u r c (r c ds,,) 
A = n, on \ r c 9A4o • 

(Note, 'the point oo' is not included in dSxo, because Qc > and A is asymp- 
totically flat, A\^ = 0.) But, in the exterior, £, ip"^ = 0, and we have the 
elliptic equation for A 

LqA = in C £ ; 

as a consequence, again by virtue of the maximum principle (Theorem ^ in 
Appendix A) now in S^q 

max A = max A = max A . 

We therefore have 

max A = max A = max A =: A{xi), for some G F . 
T F r 

Thus, (AT^.q = U£:^(,) 

maxy4 = v4(xi), for some Xi G F C A/lo (xi interior point) ; 

and, since A G C"'^(R^), in particular, A G C^{Mxq)-, it follows 

VA^ = 0. (52) 

However, this is not possible, because, on the other hand, X\ is a point of the star's 
surface Xi G F C iS, and, from the assumptions on the stellar model, dX = S 
satisfies an exterior sphere condition everywhere, i.e. dS = S U {oo} satisfies at 
each point of S (in particular, at xi) an interior sphere condition (cf. Definition in 
Appendix A). This allows us to apply the called boundary-point lemma. Theorem El 
in Appendix A, for the operator Lq in the exterior domain Sxg, with respect to the 
point xi G F C dSxo, being A{xi) = max A in S^o- And, since A ^ const, (because 
A > fic > in A/;.o, A G C^iR^), and A\^ = 0), this yields 

{v, VA)\^_^ = di^A\^_^ 7^ 0, u = outward pointing unit normal to 5 at xi ; 

contradicting (|^. Consequently, 

A<n in J. (53) 
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2nd step: A < Q < in X. 

This can be seen as follows. From inequality ()53|) and using ()48|) we also have 

Q<Qc inX, (54) 

and, combining (|53|l and (f54|l . 

A<n^ in J. (55) 

On the other hand, we have Eq. ()29|) . i.e. Lq A = —tp'^-{Q — A), satisfied everywhere, 
in particular, in the interior (in a classical sense). Let 

u{x) := A{x) — Qc Wx E I . 

Since flc is constant, we can rewrite Eq. ()29|1 . 

Lu:=LoM-^^-u = +^/'^-(r]c-^^) >0 (bydni). 

Hence, we have 

Lu>0 in J, (56) 
where, hke A, u e C\I) n C°(J), and 

M<0 in J, (57) 

by inequality We want to see m < 0. Suppose (to get a contradiction) that 

«(£) = for some x G X; then, by (|S7j), = u{x) = maXjU, x G X (interior point). 
However, by the strong maximum principle. Theorem 01 in Appendix A, for 
u cannot reach a non-negative maximum at an interior point of X, unless m is a 
constant in X. That means, in our case, u cannot vanish somewhere in X unless it 
vanishes identically in X. But u = const. = u{x) = in X, i.e. 

A = const. = fic > in X, A G everywhere , (58) 

yields, in particular, 

VA = on dI = S. (59) 

On the other hand, in the exterior, S, we have LqA = 0, with A G C'^{S) fl C°(£^), 
and, by the weak maximum principle (Theorem ^ in Appendix A) 

max A = max A , 

£ dE=SU{oo} 

but, using asymptotic fiatness {A\^ = 0) and (jKH|) . actually, 

max A = max A =: A(xi) for some Xi ^ S G dS . 
£ s 

In particular, since A ^ const., the boundary-point lemma. Theorem |21 in Ap- 
pendix A, applied to the operator Lq in the exterior domain S (where, by assump- 
tion, an interior sphere condition is satisfied in particular at xi G 5 C dS) yields a 
non-vanishing outward normal derivative 

{u,VA)l^ = d^A\^^^0, 
15 



in contradiction to (jSHl)- Therefore, m < everywhere in X, i.e. A < fic in X; 
and, hence, also on d1, because, by the weak minimum principle. Theorem ^ in 
Appendix A, applied to LqA = —ip^ ■ {Q — A) < in I (by we get miujA = 

mmgj A. Therefore, A < fic in I, or, equivalently (cf. (jlHjl ). 

A<n<n^ in X . 

3rd step: A > everywhere (i.e. the same conclusion of Proposition ^ but now using 
different hypotheses). 

We have seen in the first step A < Q in X; which yields, Lq A < in X. On the other 
hand, Lq A = in S . Accordingly, 

LqA<0 everywhere in X U £^ = . 

Applying now the strong minimum principle for generalized supersolutions. Theo- 
rem |3] in Appendix B, and using asymptotic flatness, as was argued in the proof of 
Proposition HI it follows A > everywhere. (Notice, here A ^ const., because, by 
asymptotic flatness and continuity, A = const, is equivalent to A = 0; by Eq. (jTTjl . 
also i7 = 0, and, hence, from Q — A = 0, we would have (cf. ()4fij) ) = $ = ^(0); 
but this is not possible, since requirements (i) and (ii) imply ^^-"(0) > J-'{^c) = 0.) 

Thus, A > everywhere, in particular, in the interior; using now the result of the 
second step, we finally get (jl^ . < A < f2 < f2c in X. Notice, hence, > (in X). 

4th step: A < max^ < ^Ic in £ . 

The elliptic equation holding in the exterior, 

LoA = inS, 

yields, by virtue of the weak maximum principle (Theorem in Appendix A), 

max A = max A ; 

£ de=su{oo} 

but, using asymptotic flatness {A\^ = 0), we actually have max^ujoo} A = niax^ A. 
On the other hand, we have seen A < f2 < fie in particular in dl = S, and S is 
compact. Hence, 

in < A < max A = max A < max Q < ■, 

£ S=dl S 

establishing also (fSUj). □ 

Remark 3. Analogously as argued in Remarks 1 and 2, it is possible to "reflect" Propo- 
sition 121 Accordingly, in particular, in a model for a star which is rotating differentially 
with the function JF either strictly decreasing with fic > or strictly increasing with 
Qc < 0, J-'{Qc) = 0, the angular momentum density has the same sign as Qc, and, hence, 
as the angular velocity of the fluid. Also, accordingly, the following holds 

< \A\ < \n\ < inX, 



< U| < 



maxfi 
s 



< \Qc\ in S . 
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V. GENERAL BOUNDS. ROTATIONAL ENERGY 



A. Preliminary Observation 

Let M : M" ^ M be a differentiable function, \/ : M" ^ be a vector field, and G C M" 
a domain where Gauss' theorem can be apphed. Then, due to 

dw{uV) = di{uV^ = {Vu, V)+u divV 

i 

(where (.,.) is the Euchdean scalar product, T> is the gradient operator, and divV : = 
the divergence), we get, integrating over G and applying the Gauss theorem, 

[ u diW = - [ {Vu, V)+ [ u (V, v) , (60) 

JG Jg JdG 

where u is the outer unit normal of dG (and, for simplicity in the notation, volume- and 
surface elements have been dropped). 



B. Appropriate form of the field equation 

The general (elliptic) field equation for A, Eq. may be rewritten as follows: 

diviph^^'e-^^VA) = ~f ■ (fi - A) , (61) 
where "D := {dp, d^) and 'div' are the fiat expressions in M^, and 

Jip^B 2{K-2U) 

fip, z) ^ f := pV^-^r = {e + P) ^,u_%«lu>(a-A)' ^ " ■ 

Especially we have = in the exterior E of the star (cf. ©)• 

Note, in this section (independent of Sec. IV) we go back to the field equation in 
the meridian plane coordinates, (p, z) G M.^ x M, instead of the 5-lifted one (on M^) 
(cf. Subsec. III. A). 



C. Main Observation 

Multiplying Eq. (jHH) by A, and using Eq. (jHIIl), by setting u = A, V = p^e^^e~^^VA, and 
G = ]R([x M c (actually, we consider a ball in centered at the origin of the coordinate 
system with arbitrarily large radius, B^{0) C ]R^ and take G = B^{0) n (R^ x R) c R^ 
(T — > oo), we obtain 

- f fA{n -A) = - [ p3g3Sg-4C/||p^||2 ^ f Ap^e'^e-^^'iVA, u) , 

where I C R^x R c R^ represents the (p, z)-coordinates of the interior of the star (note, 
vanishes in the exterior E). The first term on the right hand side (which converges, since 
1 1 PA 1 1 falls off rapidly enough at the space-like "infinity"'^ is obviously non-positive. And 
the second term vanishes, because of the asymptotic behavior of A at spatial infinity,^ 
and because the integrand, due to the factor p^, vanishes on the axis of rotation p = 0, 
which is the other part of 5(Rq x R) = {i? = (p^ + z^y/"^ oo} U {p = 0}. Hence, we 
have found 

jj^A{n-A)>0. (62) 
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D. Consequences 



In order to see more the linear algebra behind, we introduce now the bilinear form 



(m, v) f := J f (p, z) u{p, z) v{p, z) dpdz , m, f : / ^ M , in L (/) , 

and the induced semi-norm \\.\\f : = ((. , .)/)^/^ With this definition we can write inequal- 
ity (jni) as 

{A,n)f>\\A\\}, 

and immediately see that especially 

{n,A)f = {A,n)f>o. (63) 

Furthermore, using the Cauchy-Schwarz inequality, we have ||y4|| j||f2|| j > {A, ^)f> ll^ll/? 



and hence (since A = 
section, namely, 

i.e. 

Proposition 4 



r2 = 0) we get (in full general!) the main result of these 
0<P||^<||r]||;, (64) 




(65) 



(without any restriction concerning the rotation law, Q i-^ J^{Q) in ()14j) . in the differen- 
tially rotating case, nor the regularity and sign uniformity of These integrals can 
be regarded as some kind of "mean value" with respect to the "density" p, thus, (jU^ 
fulfilling the physical expectations.!^ 

In addition, multiplying inequality ()64|) by we find (again using the Cauchy- 

Schwarz inequality) W^Wj > \\Q\\f\\A\\f > {Q,A)f, i.e. 

{n,n~ A)f >o. (66) 

Remarkably, the integral given in ()66|) has an important physical meaning; it is, up to a 
constant factor, the called total rotational energy (see e.g. Refs. 15,16, and 18), 



T = ^j^ndj = ^ [ 2nn T^' {-gf^ dp dz = ■^{n,n-A) 



,1 °' 16 

(also denoted E-^ot or M-^ot)- Thus, (j66|) shows T > 0. Furthermore, 
16T= A)j - "^"2 / Iinii2 

by Hence, 
Proposition 5 



/ 



l^ll? 




(67) 



This generalizes the result given by Hartle (cf. Ref. 7, Sec. IV) in the limit of slow (dif- 
ferential) rotation to the general differentially rotating case. 
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VI. CONCLUSIONS 



Aiming to derive general properties of equilibrium non-singular stellar models with differ- 
ential rotation, we have established that for a wide class of rotation laws the distribution 
of angular velocity of the fluid has a sign, and then both the dragging rate (angular ve- 
locity of locally non-rotating observers) and the angular momentum density have the sign 
of the fluid angular velocity (Sec. IV). In addition, the mean value (with respect to a 
density function) of the dragging rate is shown to be less than the mean value of the fluid 
angular velocity; and this is proved in full general, without having to restrict the rotation 
law, nor the uniformity in sign of the fluid angular velocity. A further simple calculation 
of linear algebra on this inequality yields a generalization of the result given by Hartl(PI 
concerning positivity and upper bound of the total rotational energy in the limit of slow 
(differential) rotation to the general differentially rotating case (Sec. V). 
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APPENDIX A: Maximum (minimum) principles for classical sub-(super-) 
solutions 

By G we denote an open and connected set, i.e. a domain, in M", n > 2. The boundary 
is denoted by dG = G H (M" \ G). We define the differential operators 

LqU := aij{x) dijU + bi{x) diU , aij = aji , 

and Lu := Lqu + c{x) ■ u 

(where the summation convention that repeated indices indicate summation from 1 to n 
is followed), such thatP^ 

1. L (and, hence, Lq) is uniformly elliptic in G in the special form 

Q<X\y?<a,,{x)y,yj<K\y\\ Vy G \ {0}, Vx G G {\y\^ ■.= Y^y^^) ^ 

(Al) 

where A and A are constants such that < A < A < oo; 

2. all coefficients in L (and in Lq), ajj, hi (for all i and j), and c, are measurable and 
bounded functions in G, 

|ajj| < oo, < oo, |c| < oo in G (i, j G {1, . . . , n}) . (A 2) 

Theorem 1 (the weak maximum (minimum) principle for Lq (c = 0) ) 

Suppose that Lqu > {< 0) in a bounded domain G, with u G G^(G) fl G^{G). 
Then the maximum (minimum) of u is attained on the boundary, that is, 

max u = max u ( min u = min u ) . 

G dG \ G dG J 
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(A proof of that theorem can be found e.g. in Ref. 22, Theorem 3.1.) 
Definition 

For a set G C M", the boundary dG is said to satisfy an interior (exterior) sphere 
condition at a point Xi G dG iff there exists a ball B <Z G (S C \ G) with xi G dB 

Theorem 2 (the boundary-point lemma) 

Suppose that Lqu > (c = 0) m a domain G not necessarily bounded. 
Let Xi G dG be such that 

(i) u is continuous at xi, 

(ii) u{xi) > u{x) for all x E G, and 

(iii) dG satisfies an interior sphere condition at Xi. 

Then the outer normal derivative of u at Xi, if it exists, satisfies the strict inequality 

di,u{xi) > 0, 

unless u = const. = u{xi). 

(A proof of that result can be found e.g. in Ref. 23, Theorem 7, Chap. 2.) 

If c < (in Lu >Q), the same conclusion holds provided u{xi) > 0. 
(See Ref. 23, Theorem 8, Chap. 2. Also Ref. 22, Lemma 3.4.) 

Theorem 3 (the strong maximum (minimum) principle for L) 

Let Lu > {< 0) in a domain G not necessarily bounded, with u G C'^{G) fl C°(G), and 
the operator L satisfying 

c < m G (A3) 
apart from conditions ()A 1|) and ()A 2|) above. 

Then u cannot attain a non-negative maximum (non-positive minimum) at an interior 
point of G, unless u = const, in G. 

Fore = 0, i.e. L = Lq, the same conclusion holds without the requirement 'non-negative' 
('non-positive'). 

(For the proof we refer again to Ref. 22, Theorem 3.5; or Ref. 23, Theorems 5 and 6, 
Chap. 2.) 

APPENDIX B: Maximum (minimum) principle for generalized sub-(super-) 
solutions 

Consider in a domain (open and connected set) G C M" (n > 2) the differential operator 
with principal part of divergence form, defined by 

Lu = di[aij{x)djU + aj(x) u] + bi{x) diu + c{x) u , 

with = aji. Notice, an operator L of the general form Lu = aij{x)dijU + hi{x)diU + c{x)u 
may be written in divergence form provided its principal coefficients a-ij are differentiable. 
If furthermore the constants, then even with coinciding coefficients {aij = a^j, 

hi = hi, c = c) and = 0. Let us assume that 
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1. L is strictly elliptic in G, i.e. 3 a constant A > such that A < the minimum 
eigenvalue of the principal coefficient matrix [ajj(x)], 

\\y\^ <ai,{x)yiy, G M", VxgG; (B 1) 

2. Oij, Oj, 6j, and c are measurable and bounded functions in G, 

\aij\ <oo, |aj| < oo, < oo, |c| < cxd in G (i, j G {1, . . . , ra}) . (B 2) 

By definition, for a function u which is only assumed to be weakly differentiable and 
such that the functions a^djU + aiU and bidiU + cm, ? = 1, . . . , n are locally integrable 
(in particular, for u belonging to the Sobolev space Vr^'^(G)), u is said to satisfy Lu = g 
in G in a generalized (or weak) sense {g also a locally integrable function in G) if it 
satisfies 

C{u,(p;G) := / {{aijdjU + aiu)di(p — {bidiU + cu)(p}dx 
Jg 

= - [ gipdx, V^>0 ipeCl{G) 
Jg 

(where G^(G) is the set of functions in C^{G) with compact support in G). 

Notice, u is generalized suh- (super-) solution relative to a differential operator L and 
the domain G (i.e. satisfies Lu > (< 0) in G in a generalized sense) if it satisfies 
£(m, ^; G) < (> 0), > ^ G G](G). 

Theorem 4 (strong maximum (minimum) principle) 

Let u G ly^'^(G) n G°(G) satisfy Lu > (< 0) in G in a generalized sense, with 

I (c(^ - aidi'^) dx<0, > G G](G) . (B 3) 

Jg 

(equivalent to requirement ()A 3|) m i/ie classical case) and conditions ()B 1|) anc? ()B 2|) 

r/ien M cannot achieve a non-negative maximum (non-positive minimum) in the interior 
of G, unless u = const. 

(A proof of this theorem can be found in Ref. 22, Theorem 8.19.) 
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